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Compact quantum electrodynamics (CQED3) with Dirac fermionic matter provides an adequate 
framework for elucidating the universal low-energy physics of a wide variety of (2+1)D strongly 
correlated systems. Fractionalized states of matter correspond to its deconfined phases, where the 
gauge field is effectively noncompact, while conventional broken-symmetry phases are associated 
with confinement triggered by the proliferation of monopole-instantons. While much attention 
has been devoted lately to the symmetry classification of monopole operators in massless CQED3 
and related 3D conformal field theories, explicit derivations of instanton dynamics in parton gauge 
theories with fermions have been lacking. In this work, we use semiclassical methods analogous to 
those used by ’t Hooft in the solution of the U(1) problem in 4D quantum chromodynamics (QCD) to 
explicitly demonstrate the symmetry-breaking effect of instantons in CQED3 with massive fermions, 
motivated by a fermionic parton description of hard-core bosons on a lattice. By contrast with the 
massless case studied by Marston, we find that massive fermions possess Euclidean zero modes 
exponentially localized to the center of the instanton. Such Euclidean zero modes produce in turn 
an effective four-fermion interaction—known as the ’t Hooft vertex in QCD—which naturally leads 
to two possible superfluid phases for the original microscopic bosons: a conventional single-particle 


condensate or an exotic boson pair condensate without single-particle condensation. 


I. INTRODUCTION 


The parton or projective construction is one of the 
most versatile and conceptually fruitful approaches to 
a theoretical understanding of strongly correlated sys- 
tems [1]. This approach is based on rewriting microscopic 
degrees of freedom in terms of fractionalized ones that are 
charged under an emergent gauge field, and thus trans- 
form projectively under microscopic symmetries. The 
emergent gauge structure strongly constrains the low- 
energy physics, which is progressively revealed as high- 
energy degrees of freedom are integrated out. A lattice 
gauge theory with dynamical gauge fields first emerges, 
and is then replaced by a continuum gauge theory once 
lattice-scale fluctuations have been decimated. The uni- 
versal low-energy physics of the original quantum many- 
body system is then dictated by the infrared fate of this 
continuum parton gauge theory. 





Fractionalized phases of matter, such as spin liquids, 
fractionalized Fermi liquids, or fractional quantum Hall 
states, correspond to deconfined phases of parton gauge 
theories. Whether such phases exist at all for U(1) parton 
gauge theories in 2+1 dimensions—our prime focus—is a 
nontrivial question, due to the strong infrared relevance 
of the gauge coupling and the ensuing tendency to con- 
finement. Nonperturbative confinement-inducing effects 
in such gauge theories, notably monopole-instantons [2— 
4], can be suppressed by a variety of mechanisms, includ- 
ing large-flavor screening effects [5, 6], the Higgs mech- 
anism [7], and Chern-Simons topological masses [8, 9]. 
If the suppression of monopole-instantons does obtain, 
the appropriate fractionalized phase is adiabatically con- 
nected to a weakly coupled phase of the parton gauge 
theory, despite being highly nonperturbative from the 
point of view of the microscopic Hamiltonian. 


While fractionalized phases are thus perturbatively ac- 
cessible in the parton framework, conventional broken- 
symmetry phases are more difficult to describe, as non- 
perturbative confinement effects must then necessarily 
play a role. An ability to describe conventional phases 
within the framework of parton gauge theory is however 
necessary for overall consistency of the theory, as well 
as to understand the mechanism underlying confinement 
transitions between a fractionalized phase and proximate 
conventional phases. This question was studied carefully 
in recent work [10, 11] in the context of the Dirac spin 
liquid, described at low energies by U(1) quantum elec- 
trodynamics (QED3) with four flavors of two-component 
massless Dirac fermions [12-15]. Extending earlier work 
by Alicea and collaborators [16-19], Song et al. [10, 11] 
utilized the state-operator correspondence of conformal 
field theory [20] to determine the quantum numbers of 
monopole operators M for microscopic realizations of the 
Dirac spin liquid state on various lattices. The inser- 
tion of a (single) monopole operator in the Hamiltonian 
formalism corresponds to an instanton event in (2+1)D 
spacetime whereby a localized source of 27 magnetic flux 
is suddenly added to the system [2-4]. In the Hamilto- 
nian picture, a conventional phase is argued to be accom- 
panied by a monopole condensate (M) #40 which confines 
excitations with nonzero gauge charge, gives a mass to 
the emergent photon, and breaks physical symmetries if 
M transforms nontrivially under the latter. 


While these arguments are undoubtedly correct, there 
exist few explicit computations of the nonperturbative 
dynamics that would substantiate these general symme- 
try considerations. Song et al. assume a two-step sce- 
nario in which a gauge-invariant fermion mass bilinear 
first acquires an expectation value, a process described 
by an effective theory of the QED3-Gross-Neveu- Yukawa 


type [21-30] in which compactness of the gauge field is 
assumed to not play a role. After the fermionic matter 
is gapped out, instanton proliferation is further assumed 
to proceed as in the pure compact gauge theory [2-4]. 


In the presence of fermionic matter, however, gauge 
instantons may be accompanied by fermion zero modes 
(ZMs) [31], which can qualitatively affect the dynam- 
ics of instanton proliferation. Such Euclidean ZMs 
are traditionally associated with massless fermions, and 
are responsible for symmetry-breaking effects in the 
fermion sector. In (3+1)D Yang-Mills theory with 
massless fermions in the fundamental representation, ’t 
Hooft showed [32-34] that fermion ZMs on the Belavin- 
Polyakov-Schwartz-Tyupkin instanton [35] are respon- 
sible for the explicit breaking of chiral symmetry in 
the fermion sector, in a manner consistent with the 
Adler-Bell-Jackiw anomaly equation [36, 37]. Fermion 
ZMs on gauge instantons in the (2+1)D Georgi-Glashow 
model [3, 4] with massless fermions in the adjoint repre- 
sentation were shown by Affleck, Harvey, and Witten [38] 
to possibly lead to spontaneous breaking of the global 
U(1) fermion number conservation symmetry. In both 
cases, Euclidean fermion ZMs generate, via resummation 
of the semiclassical instanton gas, an effective fermionic 
interaction—the ’t Hooft vertex—that manifests the de- 
sired broken symmetry. The existence of fermion ZMs 
in the above theories is guaranteed by the Atiyah-Singer 
index theorem in (3+1)D [39] and the Callias index theo- 
rem in (2+1)D [40, 41]. The latter in particular crucially 
relies on the non-Abelian nature of the gauge field and the 
presence of a scalar Higgs field in the adjoint representa- 
tion which winds nontrivially at infinity in the instanton 
solution [3]. 

The examples above involve non-Abelian gauge fields 
and do not directly apply to our prime focus, but 
nonetheless suggest that fermion ZMs on gauge instan- 
tons may play an important role in the description of 
conventional phases and their broken symmetries in U(1) 
parton gauge theories. A natural starting point to in- 
vestigate this question is compact QED3 with massless 
Dirac fermions, relevant for the Dirac spin liquid. At 
late times and long distances, the Polyakov lattice in- 
stanton can be modeled as a Dirac monopole in 3D Eu- 
clidean space. The corresponding Euclidean massless 
Dirac equation was studied by Marston [42], but shown 
by explicit calculation to not exhibit any normalizable 
ZM bound to the instanton. Further, there appears to 
exist no generalization of the Callias index theorem to 
compact QED; [43], despite the similar infrared fate of 
the (2+1)D Georgi-Glashow model and compact QED3 
without fermions [2—4]. In the absence of explicit fermion 
ZM solutions or a general theorem guaranteeing their ex- 
istence, their relevance to the infrared dynamics of U(1) 
parton gauge theories is at best speculative. 


We emphasize here that we are interested in fermion 
ZMs bound to instantons in noncompact Euclidean 
spacetime R?, as opposed to ZMs of the Dirac Hamil- 
tonian on a 2-sphere $? surrounding a monopole inser- 


tion in the state-operator correspondence of conformal 
QED; [20]. The existence of the latter ZMs is guaran- 
teed by the Atiyah-Singer index theorem applied to the 
massless Dirac operator on the compact space S?. As 
the Marston calculation [42] indicates, however, the exis- 
tence of Hamiltonian ZMs in the latter context does not 
automatically imply the existence of Euclidean ZMs in 
noncompact spacetime. 


In this paper, we present a study of nonperturbative ef- 
fects in a U(1) parton gauge theory, in which we show by 
explicit calculation that Euclidean fermion ZMs bound 
to gauge instantons exist and lead to symmetry-breaking 
effects. The gauge theory we consider arises as the effec- 
tive continuum description of interacting lattice bosons in 
the vicinity of a multicritical point separating superfluid, 
Mott insulating, and fractional quantum Hall ground 
states [44]. While the parton description is introduced 
as a means to access the fractional quantum Hall state, 
in which a Chern-Simons term for the emergent U(1) 
gauge field leads to deconfinement, our focus here is on 
the nonperturbative gauge dynamics that obtains in the 
superfluid phase, which must result simultaneously in the 
confinement of gauge-charged excitations and the spon- 
taneous breakdown of the global U(1) boson number con- 
servation symmetry. Ref. [44] argues from general consid- 
erations that the Affleck-Harvey-Witten mechanism [38] 
should be operative and yield the desired physics, but 
does not provide an explicit derivation of the underlying 
instanton dynamics. Here we show by explicit calcula- 
tion that, by contrast with massless QED3, QED3 with 
massive Dirac fermions admits normalizable Euclidean 
ZM solutions in a U(1) instanton background. Such so- 
lutions are exponentially localized to the center of the 
instanton with a length scale inversely proportional to 
the fermion mass. Using semiclassical methods [32-34], 
we then explicitly compute the ’t Hooft vertex induced 
by those ZMs and show that it naturally leads to two pos- 
sible superfluid phases: a conventional superfluid phase 
with single-particle condensation [44], but also an exotic 
paired superfluid phase with a residual Zə symmetry. 


The rest of the paper is structured as follows. We 
briefly review Ref. [44]’s parton description of the inter- 
acting boson problem in Sec. II. In Sec. III, we formulate 
the imaginary-time partition function of the system in a 
way that makes the contribution of Polyakov instantons 
manifest, and allows us to introduce 0 parameters analo- 
gous to those of 4D Yang-Mills theory [45, 46]. In Sec. IV, 
we show that massive fermions support Euclidean zero 
modes localized on such instantons, and discuss their re- 
lationship to Hamiltonian (quasi-)zero modes in canoni- 
cal quantization. Sec. V details the calculation of the ’t 
Hooft vertex and Sec. VI explores its symmetry-breaking 
consequences. We end the paper with brief concluding 
remarks in Sec. VII; accessory technical results are col- 
lated in Appendices A and B. 


II. PARTON GAUGE THEORY 


We begin by reviewing the parton gauge theory intro- 
duced in Ref. [44]. We consider a system of charge +1 
(in appropriate units) hard-core bosons on a 2D lattice 
described by operators b(a) and bt (æ). The hard-core 
condition imposes on these operators the algebra 


[b(ar), bt (x')] = [1 — 26 (ar)b(x) 5x0’, (1) 
[b(a), b(æ')] = [b' (æ), b'(w")] = 0. (2) 


The hard-core boson then admits a parton decomposition 


b(a) = fila) falx), (3) 


where f(a) and f(x) are fermion annihilation opera- 
tors. We associate the physical boson charge with fı and 
couple this to a background gauge field A, when it is 
necessary to keep track of the physical U(1) symmetry 
associated with conservation of the boson number. The 
parton decomposition (3) also introduces a local SU(2) 
gauge redundancy f(x) > Wi;(x)f;(a), under which 
the boson operators remain invariant. In the parton ap- 
proach [1], one first ignores this gauge structure and pos- 
tulates a mean-field ansatz for the partons. Gauge fluc- 
tuations above the mean-field fermion ground state are 
then reintroduced, which ensures the parton dynamics is 
projected onto the physical boson Hilbert space. In what 
follows, we shall assume a mean-field ansatz for the par- 
tons that breaks the SU(2) redundancy down to a U(1) 
subgroup, for example via a lattice analog of the Higgs 
mechanism [7], which leaves a single gauge boson mass- 
less. Under the leftover U(1) gauge redundancy, the par- 
tons fı and fo are assigned gauge charges +1 respectively, 
so that the boson operator remains gauge invariant. 

We further consider a mean-field ansatz for the par- 
tons in which fı and f2 form independent Chern insu- 
lators with Chern numbers +1, respectively, described 
for instance by Haldane models [47] or their analog on 
the lattice of interest. In the vicinity of Chern-number- 
changing transitions in the parton bandstructure, this 
theory is described in the continuum limit by a 3D Eu- 
clidean Lagrangian, 








L= DD [pial = iA + m)Yia + Pal a m)2e| , (4) 





where A, is the background field that tracks the physi- 
cal U (1) symmetry, and {W14, Y2+} are two-component 
Dirac fermions obtained in a linearization of the partons 
{ fı, fo} at the two Dirac points K+ that generically ap- 
pear in the parton bandstructure [48]. Importantly, the 
fermion masses for 1+ and 24 are opposite in sign, 
since the Chern numbers are opposite in sign for fı and 
fo. 

For the above mean-field parton ansatz to correspond 
in fact to a physical state of bosons, we must reintro- 
duce gauge fluctuations. To study the effect of those 
fluctuations, the lattice fermions (i.e., the partons) are 























minimally coupled to an emergent U(1) gauge field 
Gy. For example, the parton fı of gauge charge +1 
minimally couples to the gauge field on link (a#,7) as 
Fates exp(—ida,i) fix+i, where tæ, is a hopping integral, 
z is a lattice site, and 7 is a lattice vector. The invariance 
of such a term under 2nz shifts of az; can be viewed as 
a gauge redundancy or as a true local symmetry. These 
two perspectives will be discussed in Sec. III. In either 
case, at low energies, the renormalization group endows 
the emergent field a with dynamics that preserves this pe- 
riodicity, implying an effective gauge field Hamiltonian of 
the form 

















y= 5d + KY cos ), (5) 











where e is the electric field on link | satisfying [a;, ey] = 
tou, and fo is the flux (lattice curl) of a through the 
plaquette O (we shall henceforth assume a square lattice 
for simplicity). The physical Hilbert space of the gauge 
theory (with fermions) is the gauge invariant subspace 
specified by a Gauss constraint. Weak fluctuations of 
a, correspond to the K > 1 limit, in which H, is en- 
ergetically appeased by fg = 2ang, where ng € Z is a 
plaquette-dependent integer. Expanding about any one 
of these minima leads to the usual Maxwell theory with 
a massless photon. However, it is well known that tun- 
neling events fg > fo+27Q, where QEZ, on a plaquette 
cannot be ignored, for these give the photon a mass expo- 
nentially small in the gauge coupling K. These tunneling 
events, corresponding to 27Q flux insertions on a plaque- 
tte, feature as instantons (Dirac monopoles of charge Q) 
with finite action in the 3D Euclidean theory [2-4]. 

In a naive continuum limit, the effective parton La- 
grangian with gauge fluctuations is 










































































L= X [Pra(d — ih- id + mre 





+ÜzalÖ + if = maa] + pal? (6) 


where e is the renormalized gauge coupling (some func- 
tion of the lattice coupling K). However, a finite UV 
regulator (lattice constant) and the fact that the lattice 
theory is invariant under a— a42nz imply that the effects 
of instantons must be accounted for in this continuum 
limit. This theory is termed compact QED; (CQEDs). 
We note however that by contrast with the CQED3 the- 
ory of the Dirac spin liquid, which also has four flavors of 
two-component Dirac fermions, the fermions in our case 
(i) are massive, and (ii) do not all have the same sign of 
the gauge charge. 


III. 0 PARAMETERS AND INSTANTONS 


In this section, we use canonical quantization to derive 
a path integral representation of the partition function 


of the pure gauge theory without matter, which makes 
the contribution of instantons explicit and allows us to 
introduce 0 parameters [49, 50] analogous to those of 4D 
Yang-Mills theory [45, 46]. This sets the stage for our 
computation of the ’t Hooft vertex using path integral 
methods in Sec. V, after explicit fermion ZM solutions 
in the background of a single instanton are obtained in 
Sec. IV. 

We begin with the pure gauge theory, described by 
the Hamiltonian (5), which we shall consider in the ab- 
sence of background charges. This means that the Gauss 
constraint on every site is (dive), |Y} = 0 on all phys- 
ical states |W). As stated previously, the invariance of 
H, under 27Q translations of the gauge flux on a pla- 
quette can be viewed as either a true local symmetry, or 
as a gauge redundancy due to rotor-valued link variables 
awi E U(1) = R/2nZ. The former view will be called 
minimal compactness, and the latter forced compactness. 
In what follows, we shall mostly be concerned with the 
“magnetic limit” K >> 1, in which gauge fluctuations are 
weak. 


A. Minimal compactness 


In the minimal compactness picture, the gauge field 
axi ER. A general state in the Hilbert space is given by 
a wavefunctional 


Vlasi] = ({aea}|¥), (7) 


where {aw i} denotes the collection of a on all links, and 
\{az,:}) forms a basis. The electric fields ey, generate 
translations of these wavefunctionals. On a single link, 


e= a) = |a +a), (8) 


which means 


A gauge transformation exp[—id(dive),] on a site æ is a 
translation that leaves all plaquette fluxes fg invariant. 
Since the magnetic term is a periodic function of fo, Hg 
is not only invariant under these gauge transformations, 
but also under a discrete group of local flux translations 
fo- fo+27Q, QEZ on a plaquette. This group is gen- 
erated by monopole operators Mg(#), where & denotes 
a plaquette (or equivalently, a site on the dual lattice), 
and 












































M} (@) feMa(#) = fa + 27Q. (10) 


This translation is also generated by electric fields, but 
one must use an infinite string of fields [51], since only 
the flux in plaquette z must be changed. One possibility 


is to consider an infinite product of all horizontal links 
below z, and non-uniquely define 


0 
DD CEx+pĉz2,ĉ13 (11) 


p=—0o 


Ma(#) = en, Oz 


where ĉ1,ĉ2 are unit vectors in the positive horizontal 
and vertical directions, respectively. 

The minimally compact theory has similarities with 
the Bloch problem of electrons in a crystal lattice, in 
which a discrete translation by a lattice constant is a 
physical symmetry as opposed to a gauge redundancy. 
In the Bloch problem, there occur instantons that tunnel 
between the minima of the crystal potential, and the true 
ground state is a superposition of all local minima. In 
minimally compact CQED3, the analogs are monopole- 
instantons that tunnel between physically distinct min- 
ima fz=2n7 to f,=2m7 on a given plaquette %. Since 
[H,, Mg(#)|=0 on every plaquette, the physical eigen- 
states of H, fall in representations of these symmetries. 
The irreducible representations of this Abelian group are 
all one-dimensional, and are simply phase factors (Bloch 
theorem). The eigenstates of H, must thus obey, 

Va, Ma(#) |) — aks |Wn,o) ’ z€ [0, 27), (12) 
where @z is an analog of crystal momentum, and n is 
a collective index denoting all the other quantum num- 
bers necessary to specify the state. The corresponding 
eigenenergies will be continuous functions of 0z, as in 
conventional band theory. 

For example, a single square plaquette in the “electric 
limit” K — 0 (the analog of the “empty lattice approxi- 
mation” in the Bloch problem) is governed by the Hamil- 
tonian 
1 
H, = z ezi 


wi 


(13) 


There is a single site & on the dual lattice, and we thus 
drop the site index. Eigenstates of H, are eigenstates of 
all four electric fields bordering the plaquette, but sub- 
ject to the Bloch condition (12) and the Gauss constraint 
(dive), |¥)=0. For this single-plaquette system, Eq. (11) 
implies og = €,¢, and thus Mg(#) = exp(t27Qez, 3, ). 
The Bloch condition (12) is then 


eTRes, i |Vno) = e28 |W.9) ‘ (14) 


This implies physical eigenstates of ex, (and of H, in 
the limit K — 0) satisfying the Bloch condition are re- 
stricted to those with eigenvalues 


neZ. (15) 


Ex, SN + On’ 
The electric fields on the other links can be found using 
the Gauss constraint. The physical states are loops of 
electric flux circling the plaquette, with an integer level 


spacing. Substitution of these values into the Hamilto- 
nian (13) gives the bandstructure 


E,(0) =2 (n+ ae (16) 


In the minimal compactness picture, {0z} are a set of 
quantum numbers specifying states in the same Hilbert 
space. 

Finally, we observe that [H, Mg(#)|=0 for the full lat- 
tice Hamiltonian H with gauge fields and fermions, since 
the gauged fermion hopping term discussed in Sec. II is 
invariant under local shifts of the link field ag, by arbi- 
trary integer multiples of 27, including those produced 
by conjugation with the monopole operator (11). Thus 
the Bloch condition (12) applies to eigenstates of H as 
well. As with the Bloch theorem in solid-state physics, 
an eigenfunctional of H satisfying this Bloch condition 
can be written as the product of a “plane wave” and a 
periodic function, 


Wn olaa] = er Xa lafa Pn olaz,il, (17) 


where ©, 9 is invariant under 27Q flux translations on 
a plaquette, i.e., Mg(%)® =, and we have suppressed 
the dependence of ® on fermionic coordinates, which does 
not play a role in this analysis. As in band theory, we 
can reduce the solution of the Schrödinger equation over 
ax i E€ R to that over a single “unit cell” az € [0, 27) 
by either solving the original equation over that domain 
with the twisted periodic boundary conditions (12), or by 
deriving an equation for the periodic part ®. Defining the 
unitary 


Up = e% Xa lafa, (18) 


we see that ®,, 9 obeys the modified Schrödinger equation 
Ho®n o =En(0)Pn,o, where 


Ho = U} HU 


1 1 2 
P2 feast zesta, +K X (1-cos fa). 
(19) 


II 


Since 0% only enters Ho through its spatial lattice deriva- 

tive, a uniform parameter 0z = @ has no effect in the 

bulk [49, 50], but will affect energetics in a system with 

boundary as the single square plaquette considered here. 
The partition function for a fixed set {0%} is [52] 


Zo = tr e PHe — 5 e Le Sata Zo, (20) 
{Qa}EZ 


where the second equality is a Fourier decomposition, 
since Zo is periodic in all the 0z. The Fourier coefficients 
are given by 


2m 
Zo= | D{Oz}e~* Xa 22% Ze, (21) 
0 


5 


where we define D{6z} = [], %2. These Fourier coef- 
ficients can be interpreted as a partition function of the 
original Hamiltonian H with monopole insertions as fol- 
lows. Let M(Qz) be a product of monopole operators 
that inserts flux across the system in a manner deter- 
mined uniquely by the configuration function Qs. Using 
the completeness of flux (f) eigenstates in the gauge- 
invariant subspace, we obtain: 


treP™M(Qg) = 1 D{ fa} {fa} e PE M(Qa) fa} 


z [ DU fo} Ufot?" fo + 27 Qg}) 
(22) 


df; 
where D{ fg} = Tl, Z. 
eigenstates of H using 


Inserting a complete set of 


I= f Pie} 37 Wao) (nol, (23) 


we find that 


rer" 3) = 7 z 5 
tre~>¥ M(Qa) i D{0z} L Pm E 


x ({fg tle 7 (Eno) (Vn olfa +27Q3}) 


2r 
=} D{0a} Sy EPPO ei Ea Qala 
0 
{na} 


x 1 Df fy} T3 olflYn olf] 


27 
= D{0z} e~ Le Qaba ty e—PHe 
0 
= ZQ. (24) 


In the second line, we have used the Bloch condition sat- 
isfied by the gauge-invariant wavefunctional VW, [f] as 
seen in Eq. (17), and the third line follows from its nor- 
malization to unity. Therefore, the partition function for 
a fixed set {0z} can be written as 

Zo = 5 e La Qabe tre PH M(Qz). (25) 
{Qa }EZ 


Each term of this series can be written as a path integral 
in a fixed monopole configuration background. The Q- 
dependent exponential prefactor can be absorbed into the 
trace by explicitly including a 6-term i f 0(a)sdf (a) in the 
action, where f(x) is now the Euclidean electromagnetic 
2-form and x denotes the Hodge star. 


B. Forced compactness 


In the forced compactness picture, the gauge field 
axi € U(1) = R/27Z is a rotor-valued variable. The 


canonically conjugate electric fields then have a spectrum 
valued in Z. In this perspective, the various minima 
fz = 2nz7 for ng € Z are identified as the same state, 
and a flux translation fz —> f2+27Q becomes a gauge 
redundancy. 

In this perspective, the problem is akin to that of a 
quantum particle on a ring, where a translation by a 
length equal to the circumference is a gauge transforma- 
tion. If the ring is suspended in a gravitational potential, 
then the unique classical ground state that minimizes 
the potential is at the bottom of the ring. However, in 
the quantum problem, there are tunneling events (instan- 
tons) that correspond to the particle winding around the 
ring an integer number of times, which involves overcom- 
ing a potential barrier. 

The exact analog in CQED3 in the forced compactness 
picture are monopole-instantons that cause fz > fzgt27Q 
on a plaquette. The monopole operator defined by 
Eqs. (10)-(11) is thus a gauge transformation (a “do- 
nothing” operator) that connects different labels for the 
same physical state. These are called large gauge trans- 
formations, terminology inspired by analogous concepts 
in 4D Yang-Mills theory [45, 46]. Large gauge transfor- 
mations are distinguished from the usual small ones in 
that the former crucially utilize the multi-valuedness of 
the gauge function. A small gauge transformation is one 
of the form [[„ explida(dive)s], where {ġa} are single- 
valued gauge functions, i.e., all of them lie in a single 
branch of R/27Z, for example [0,27). One example of 
a large gauge transformation is Mpo exp[20(dive)2+p2, | 
which, since 0 ~ 27 in R/27Z, can be written as 
Ilpo exP(i27€24p%1,22), Which is the same operator as 
(11), but with a displaced Dirac string. 

Physical states are required to be invariant under all 
gauge transformations, small and large. This imposes 
0z =0 for all plaquettes z in the Bloch condition (12). 
For the particle on a ring, a background flux can be 
threaded through the ring, which changes the Hamilto- 
nian and the Hilbert space of the problem, as we are deal- 
ing with a physically different system. The background 
flux can be unitarily removed from the Hamiltonian, at 
the expense of twisting the boundary conditions on wave- 
functions, which in winding around the ring, will then 
gain an Aharonov-Bohm phase. Similarly, in CQEDs, 
one can introduce a theta term, which changes the Hamil- 
tonian from Eq. (5) to Eq. (19). There is a macroscopic 
number of such possible theta terms, corresponding to 
a choice {0z}. Again, one can remove the theta terms 
from the Hamiltonian, but at the expense of introducing 
twisted boundary conditions under large gauge transfor- 
mations, as in Eq. (12). The key difference with the min- 
imal compactness picture is that a given set {0z} labels 
the entire Hilbert space of the theory, sometimes called 
a given theta universe. States with different {0z} belong 
in different Hilbert spaces; conversely, states in the same 
Hilbert space have the same {0z}. 

The expression (25) for a partition function with fixed 
{0z} remains valid in the forced compactness perspective. 


In fact, it is the full (i.e., unrestricted) partition function 
here since the entire Hilbert space is characterised by the 
fixed set of parameters {0z}. 


IV. ZERO MODES OF MASSIVE FERMIONS IN 
INSTANTON BACKGROUNDS 


Having discussed instantons in the pure gauge the- 
ory, we now include fermions. As mentioned previously, 
the presence of fermionic ZMs in instanton backgrounds 
is typically associated with massless fermions in non- 
Abelian gauge fields [32, 33, 38, 40, 41]. Marston [42] con- 
sidered massless Dirac fermions in an Abelian instanton 
background in (2+1)D and found no fermion ZMs bound 
to the instanton. Motivated by the U(1) parton gauge 
theory (6), we consider here massive Dirac fermions in 
the same background and show by explicit construction 
that fermion ZMs now exist. This result is in accordance 
with the existence of zero-energy bound states for rela- 
tivistic fermions in a (soliton) monopole background in 
(3+1)D (see Ref. [53] and references therein). In the soli- 
ton version of the problem, the fermion ZM is found by a 
self-adjoint extension of the Dirac Hamiltonian. Such a 
technique is inapplicable for the instanton version of the 
problem as the Euclidean Dirac operator D appearing in 
the action [see Eq. (26)] is not Hermitian, nor is there any 
requirement for it to be. Rather, D must obey reflection 
positivity (see Appendix B). Therefore, the calculation 
of the ZM solution must be done anew in the context of 
the instanton problem. 

As Marston himself notes, the Callias index theorem 
for odd-dimensional noncompact manifolds provides the 
number of fermion ZMs in the case of massless fermions 
in the background of non-Abelian instantons [40, 41]. 
This index theorem crucially relies on (i) the existence 
of a Higgs field, and (ii) on relating the index of the 
Dirac operator D to (dim ker DİD — dim ker DD"), both 
of which fail to hold in the current setting of massive 
fermions in Abelian instanton backgrounds. The reason 
for the failure of (ii) might seem surprising, and is dis- 
cussed in Appendix B. Despite the absence of a rigorous 
index theorem for the current problem, and as we discuss 
below, the ZMs we find by explicit solution can be given 
a topological interpretation by analogy with Hamiltonian 
(quasi-)zero modes associated with the Atiyah-Singer in- 
dex theorem. 





A. Setting up in spherical coordinates 


A 3D Dirac fermion w of charge e (with sign) and mass 


m, in an instanton background a% ) of topological charge 


g, is specified by the Euclidean Lagrangian 
L= 4 — ieg + m)p = YDY. (26) 


The instanton is assumed to sit at the origin, with its 
Dirac-monopole vector potential defined à la Wu and 


Yang [54]. Working in spherical coordinates (r, 0, p), we 


have al) = al?) =0 and 


al?) _ 


where charts in the northern and southern hemispheres, 
Ry and Rgs, are defined by Ry : 0 € [0,7/2+6) and 
Rs : 6€ (1/2—6,7], and a choice of 6 € [0,7 /2) defines 
the chart overlap region Rug : 0 € (7/2—6,7/2+6). 
The Dirac matrices are Pauli matrices (Ox, Cy, oz), with 
2 being the Euclideanized time direction. The ZM wo of 
the Dirac operator D solves 





(cosð — 1), re Ry, 


cos +1), r€ Rs, 


g 
rsin 


(27) 





rana 


iDuo = (-i9 — ed 





im)o = 0. (28) 


We will treat this formally as a quantum mechanics prob- 
lem in three spatial dimensions, regarding —i0; as a 
canonical momentum operator p;. The ZM equation is 
then 


(o;p; — egoja? 


where the mechanical momentum 7 


p—ea\), Defining 
gf =o -7, we use the fact that (o - r= 


o2=1 to write 


ft =(o-t)(o-f)(o- 7) 


=(o-f)[f-m+io- (ê x r). (30) 


The canonical angular momentum that is conserved in a 
monopole field is (see Appendix A) 


L=rxnrn -— egf, (31) 


where eg €Z/2 by the Dirac quantization condition. This 
can be used to rewrite Eq. (30) as 


t= (o-i) (trti got), (32) 


Since L generates spatial rotations, which leave r = |r| 
invariant, [r, L]=0 and the placement of r does not mat- 
ter in the above equation. Using m =p—ea') and since 
f- al?) =0 for the Wu-Yang potential, 


ft =(o-f) (t p+io r4 “fo-2) 
r r 


= ~i(e -#) (a -io-r- 20.8), (33) 


The Dirac operator is thus 








ae Za) Hm. (34) 


B. Zero modes of the Dirac operator 


Firstly, we note that [J,D]=[J,i7#]=0, since f=o-7 
is a dot product that remains invariant under simultane- 
ous rotations of ø and m generated by the total angular 
momentum J = L + io. As a set of commuting observ- 
ables, we take 


[J?, D] = [Jz, D] = 0. (35) 


This means the angular part of eigenspinors of D are 
the monopole spinor harmonics DV jmi (8, p) (see Ap- 
pendix A), which informs the eigenspinor ansatz: 


Wim; = Aga, (VAN? (0, P) + Bjms VE, (8,9); 


(36) 
for a given (eg,j,m,). It is necessary to superpose both 
values of the orbital angular momentum, L+ = j +1/2, 











that give rise to a given total j as L is not a good quantum 
number. 

For a monopole of the lowest magnetic charge g = 
+1/2e, the total angular momentum can assume j=0, 1. 
We will focus on the lowest spherical wave j =0, for which 
the orbital angular momentum L_ = —1/2 is excluded. 
To reduce notational clutter, we suppress (eg, j, mj) la- 
bels everywhere except on the monopole spinor harmon- 
ics, and write the ZM ansatz as: 














ol, 0, p) = Ax (r)Y 72,0009 9): (37) 





The zero index on wo indicates that it is a ZM. The 
indices on the radial part A4 (r) denote the value of eg 


+1/2. Finally, Yizo, stands for the (L = 1/2,eg = 
+1/2, j = 0,m; = 0) monopole spinor harmonic. Since 
this ansatz is coincidentally also an eigenstate of L?, we 


rewrite the Dirac operator in Eq. (34) as: 


I H 

















D= (a-i) |ð, L(a r >) “oy «| +, (38) 


The action of D on the ZM ansatz (37) is then (for eg= 


+1/2), 





is 3 
(a -f) Piano +54 w 

















2r = 





eo Alo fyi? 


or +1/2,0,0 





+ MALY i999 =0. (39) 




















Using Eq. (A22) of Appendix A with eg=+1/2 yields: 











(o. A i E E ot: (40) 


and therefore: 














BAIA G + m) Ai 20, (41) 


which have the obvious exponential solutions. Therefore, 
the ZMs for eg=+1/2 are 


2m -mr 
y(r 0 p) = He Ya 0,0(8, 9); 


y —2 
Memryi/2 
r 





Yo (i; 0, p) = —1/2,0,0 (9, p), (42) 
which are normalizable for m>0 and m <0, respectively 
(recall that there sits an r? in the Jacobian for spherical 
integrations). The divergence at r=0 is superficial. The 
origin is excluded in the problem due to the instanton 
there, mathematically implemented by working in spher- 
ical coordinates. Alternatively, we know that the instan- 
ton has a finite core due to the underlying lattice, and 
so the derived form of the ZM is valid only at large dis- 
tances compared to the lattice constant. This has been 
discussed further in the soliton version of the problem by 
Yamagishi [53]. 


C. Zero modes of the adjoint Dirac operator 


Under an integration by parts, the Lagrangian changes 
from YDy to (Dypt) tp, so the ZMs of the adjoint Dirac 
operator are also important. The adjoint Dirac operator 
is 


Dİ = -ð+ ied +m, (43) 


where the adjoint is defined with respect to the standard 
inner product on L?(R*). In fact, the correct domains 
of D and DÝ are subsets of L?(R*), as discussed in Ap- 
pendix B. The ZM equation for D? is 


iDi ho = (-id — ed + im) = 0, (44) 


which is the same as that of D, but with the sign of m 
reversed. This implies the ZMs of DÝ for eg =+1/2 are 








T vV—2m mr 

Bo (r, 0, p) =e" Via 0,009), 

p V 2m —mr 

Wo (r, 0, p) = aan C am (9, p), (45) 


which are normalizable for m <0 and m >Q, respectively. 
We note that D and DÌ cannot both possess normalizable 
ZMs simultaneously. 


D. Hamiltonian picture and quasi-zero modes 


Although the Callias index theorem does not directly 
apply to the problem considered here, the topological 
nature of the fermion ZMs explicitly found in the pre- 
vious subsections can be understood by considering an 
approximate treatment of the same problem but in the 
Hamiltonian picture, following the approach of Refs. [16- 
19]. In this approach, a single instanton is modeled as 


a static, infinitely thin 27 solenoidal flux to which the 
fermions respond: V x a = 2rô(r), where r= (x,y) now 
denotes the spatial coordinate. Consider first fermions 
of gauge charge +1 and zero mass. As is well known, 
for each fermion flavor, the corresponding 2D massless 
Dirac Hamiltonian possesses a single quasi-normalizable 
“chiral” zero-energy eigenstate Wo = (u,0)? with u ~ 
f(a+iy)/r, which can be understood as a manifestation 
of the Atiyah-Singer index theorem [55]. 


In the presence of a nonzero mass m > 0, a “zero-mode” 
solution persists, again of the form Yọ = (u,0)? with 
u~ f(a+iy)/r, but its energy E is shifted from zero to 
E =m. For fermions of gauge charge —1 and negative 
mass —m, there is also a single such quasi-zero mode 
per fermion flavor, now with wavefunction Yọ = (0,v)" 
and v ~ g(a—iy)/r, but again energy E =m. For a 
—2r flux background corresponding to an anti-instanton, 
the situation is reversed: fermions of gauge charge +1 
and mass m possess a quasi-zero mode w= (0, v)? with 
v ~ g(x—iy)/r, and fermions of gauge charge —1 and 
mass —m possess a quasi-zero mode Wo = (u,0)" with 
u~ f(a+iy)/r, both with energy E=—m. 

The Hamiltonian quasi-zero modes discussed above are 
similar to those appearing in the “zero” mode dressing 
of monopole operators at the quantum critical point be- 
tween a Dirac spin liquid and an antiferromagnet [24, 56, 
57]. In the latter context, a spin-Hall mass mo, appears 
spontaneously in the saddle-point free energy of the asso- 
ciated conformal field theory quantized on a sphere sur- 
rounding a monopole insertion, following the approach of 
Ref. [20] to calculate the scaling dimension of monopole 
operators at conformal fixed points. This spin-Hall mass 
gives a mass of opposite sign to fermion flavors y and 
p, of opposite spin, but the gauge charge is the same 
for both flavors. In our case, the two parton flavors %1 
and wz in Sec. II play the role of yp and y, and the 
“spin-Hall” mass comes from the parton Chern numbers 
appropriate to a superfluid state [44]. Furthermore, the 
gauge charge is opposite for both flavors on account of 
the parton decomposition (3). Nonetheless, in both cases 
a single normalizable quasi-zero mode with energy +m 
appears for each Dirac fermion flavor, as expected from 
the Atiyah-Singer index theorem. 





Finally, the counting of instanton zero modes in 
Sec. IV B and IV C is consistent with that of the Hamil- 
tonian quasi-zero modes just discussed, if both D and 
Dt are considered. For an instanton (g > 0), and for 
fermions of positive gauge charge (eg = 1/2) and mass 
m>0, the Euclidean Dirac operator D has a single nor- 
malizable zero mode og xe" /r, where r now denotes 
Euclidean spacetime distance from the center of the in- 
stanton [see Eq. (42)]. For fermions with negative gauge 
charge (eg =—1/2) and mass m<0, D likewise possesses 
a single normalizable zero mode, Wg x e™"/r. The ad- 
joint Dirac operator Dt has no zero modes in this case. 
For an anti-instanton (g < 0), the situation is reversed, 
as in the Hamiltonian picture: D? now has normalizable 
zero modes, but D has none. For e> Qand m>0, Dİ has 


a single zero mode Wo xe” /r; for e<0 and m<0, the 
Dt zero mode is of xe” /r [Eq. (45)]. The fact that in- 
stantons (anti-instantons) are associated with zero modes 
of D (DÌ) is further discussed in Sec. V B and V C and has 
important consequences for instanton-induced symmetry 
breaking. 


V. THE ’T HOOFT VERTEX 


Using the explicit ZM solutions for massive 3D Dirac 
fermions found in Sec. IV, we now show that Abelian 
instantons in CQEDs3 induce symmetry-breaking inter- 
actions for such massive fermions. The calculation here 
is analogous to ’t Hooft’s groundbreaking solution of the 
U(1) problem in 4D quantum chromodynamics [32, 33], 
which is well summarized in a review article [34] by the 
same author. In short, fermion ZMs cause transition am- 
plitudes with nonzero instanton charge Q £0 to vanish 
when evaluated between states of equal fermion num- 
ber. Instead, nonvanishing amplitudes occur between 
states of different fermion number: fermionic field inser- 
tions appearing in such amplitudes “soak up” the fermion 
ZMs appearing in Grassmann integration. Resumming 
these insertions in the Coulomb gas formalism produces 
a fermionic interaction term known as the ’t Hooft ver- 
tex, whose symmetry is lower than that of the classical 
Lagrangian. 


A. Partitioning the partition function into 
instanton sectors 


To make the calculations less tedious, we consider only 
two species of 3D Dirac fermions %1, Y2 with opposite 
gauge charges e; = —e2g =e and zero net Chern number, 
so that their masses satisfy mı = —m2z =m. This cor- 
responds to effectively ignoring the valley (+) subindex 
in the original Lagrangian (6), which can be easily re- 
stored at the end of the calculation (see Sec. VI). The 
Lagrangian of interest is then 


L= (9 — id +m) + p(l + id — my 


1 ; 
+ga fu + 40(x)ewwrOufra, (46) 








where an explicit 0 term has been added to keep the dis- 
cussion general (see Sec. III). The fermion part of the 
action will be denoted as Sp[a,,] when there is need to 
refer to it separately. The presence of a lattice regula- 
tor permits monopole-instantons in this theory to have 
finite action. The qualitative effects of those instantons 
on fermions are what we wish to understand. To for- 
mulate this problem in terms of a path integral, we use 
Eq. (25) to write the partition function in a fixed 0 uni- 
verse as 


Z=5 | Data) Dange, (47) 


QEZ 


where [Da,]o indicates a restriction of the integration 
over a, to configurations with total instanton charge 
Q/2e. The sum over instanton configurations in Eq. (25) 
has been replaced by a sum solely over total charge Q, 
with integrations over instanton collective coordinates 
subsumed in the measure [Da,]g. We will “integrate 
out” the instantons in the Q 4 0 sectors and write a 
theory of fermions. A direct coupling between photons— 
quantum fluctuations of the gauge field about the classi- 
cal instanton background—and fermions remains in this 
final theory, but can be neglected in the computation of 
the ’t Hooft vertex, which is a semiclassical effect [32-34]. 

Insofar as the pure gauge theory is concerned, the effect 
of an instanton of charge Q/2e is captured by an insertion 


ee? J Bx e2t (2), (48) 


in the path integral [2-4]. Here, t?Q?/2c? is the action of 
a charge-Q/2e instanton, o(a) is a compact scalar called 
the dual photon—the continuum limit of the lattice vari- 
able og in Eq. (11)—and the integration is over the in- 
stanton position x. Likewise, the factor exp(i27Qc) is 
just the spacetime representation of the monopole oper- 
ator Mg [20], whose lattice form was given in Eq. (11). 
An attempt to dualize this theory in the same vein as 
the classic Polyakov duality between the compact U(1) 
gauge theory and the sine-Gordon theory [2-4] leads to 
the path integral 


3 ; lore) 

z -F S Pao)? 5 a 

Z= | Doe 2 NI 
N=0 


2 


N co 
X [I [ez NO eT EOR eiQrl2no ln) +O] 
k=1 Q 


x I D(a, paje” 5r", (49) 


There are a number of results used in writing Eq. (49), 
especially with fermions present. The gauge field in each 
Q-sector in Eq. (47) has been formally decomposed as 


ap = a3 W + aQ, (50) 


where a W is the photon (analog of “spin wave” in the 
2D XY/sine-Gordon duality [58]) part whose coupling to 
fermions is neglected at the semiclassical level, and ag de- 
scribes an instanton configuration of total charge Q/2e. 
The photon part has been dualized to the Gaussian ac- 
tion for the compact scalar ø. In the above decomposi- 
tion, the photon part can be thought of as finite-action 
fluctuations (of arbitrary size) around the fixed instan- 
ton solution. In addition, the sum over charges Q, and 
integration over positions implicit in the measure [Day] 
in Eq. (47) have been rewritten as sums over the number 
N of instantons, their charges Qk, and their locations z,. 

An important assumption used here is the dilute gas 
approximation, which gains new significance in the pres- 
ence of fermions for the following reason. Consider, for 


instance, N = 2 with charges Q/2e and Q2/2e. It has 
been assumed [32-34] that the gauge field for such a con- 
figuration can be written as 


aQ1T92 (x; z1, z2) ~ a (x; z1) + a2? (x; z2), (51) 


with |z1 — 22| 1. This might seem plausible given the as- 
sumption of a dilute instanton gas, but the consequence is 
severe, for it implies that the Dirac action also separates, 


Splagi*@] ~ Spla?] + Sr lag’). (52) 


This allows the fermion path integral to be written inside 
the N and Q sums, as in Eq. (49). If Qı =-Qə2 so that 
the total instanton charge is zero, one might expect that 
there exist no normalizable fermion ZMs. However, the 
decomposition of the action in the above form (which fil- 
ters into a decomposition of the Dirac operator) clearly 
allows ZMs. We shall nevertheless assume such a de- 
composition; the error in the resulting partition function 
turns out be of O(A?), where 


Aer, (53) 


is the action for an instanton of lowest charge, and we 
consider the semiclassical limit A« 1. Moreover, even if 
there are no strict (topologically protected) ZMs in this 
case, there will likely be eigenmodes of the Dirac operator 
lying arbitrarily close to zero, with splitting proportional 
to exp(—m|z1 — 22|), since the ZM wave functions (42) for 
an isolated instanton decay exponentially away from the 
center of the instanton. 

In what follows, we shall only consider the effects of 
a dilute gas of instantons of elementary charge, thereby 
restricting the sum over Q; in Eq. (49) to +1. The con- 
tributions of instantons of higher topological charge are 
suppressed by increasing powers of A. 





B. Q=1 instanton sector 


In the Q=1 sector, the fermion part of the partition 
function is 


Zrlaj] = / D(a, hae SF, (54) 


where at describes a single charge 1/2e instanton located 
at z}. We will show that this path integral is precisely 
zero. Before the zero is revealed, the fermion functional 
measure needs definition. For technical reasons discussed 
in Appendix B, the standard means of definition using 
the eigenfunctions of DD and DD? does not work, since 
these only span subspaces of L? (R) over which DID and 
DD! are self-adjoint (not merely Hermitian), and it so 
happens that the ZM lies outside these subspaces. In- 
stead, we shall proceed along physical lines; the effects of 
ZMs are what we are interested in. The results of Sec. IV 
indicate that the Euclidean Dirac operators, 


p tM, 


DSO, + iqta, — m, (55) 


Dı = y"0, — iy"a 
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each possess a normalizable ZM. We shall thus use a 
mode expansion of the Fermi fields as 


yı (z) = wo(@ — z+) + py) ui(@ — 24)ni, 


7 


hlz) = vole = z4)xo + JD vile- z4)xi (56) 


a 


where ni, Xi are single-component Grassmann variables. 
ug and vo are the respective ZMs of Dı and D2, calcu- 
lated in Sec. IV B, and the primed sum denotes non-ZM 
contributions. As discussed in Appendix B, we can ei- 
ther assume that there exists some self-adjoint operator 
whose domain includes the ZM, or we can use the eigen- 
functions of a self-adjoint Dİ Da (which excludes the ZM) 
to account for non-ZM contributions and add the ZM by 
hand. Either way, the ZM contribution has to be ac- 
counted for on physical grounds. 

Since Di and DÌ are deprived of normalizable ZMs in 
the Q = +1 sector, we write mode expansions of 1 € 
Dom(D!) and p2 € Dom(DÅ) as 

— / 

Yi (x) = 5 ti(z4 — £) Tlz4 — x) = ù] (x — z4), 
(3 

= / 

palz) = D> (zy -2)Xi, T(z4 — 2) =ò] (z — z4), 


i 








(57) 


where the transpose acts on spin and spatial indices, no 
matter where the T is placed. The functional measure for 
fermions can now be defined as 


= = i 
D(b1,1, P2, 2) = dno dxo | | d7: dn: dxi dx: 


= dmo dxo D'(j,n)D'(X,x). (58) 


Since the ZMs nņo, Xo do not appear in the action, the 
path integral for Zp [Az] is killed by the measure. There- 
fore, only the sector with zero instanton charge con- 
tributes to the full partition function (Z) of the theory. 

However, sectors with nonzero instanton charge will 
contribute to correlation functions that can “soak up” the 
ZMs. For instance, (only) the Q=1 sector contributes to 
the correlation function 


- Z( Puy («)a(22)) 
x = l dno dxo D' (ñ, n, X, x) > uT (£1 — z4 )ni 


x 5 vj(z2 — 24) xje7 Sr laa 
j 


_ a z at 
= ug (ay — 24)00(02- 24) | D'O 2de Sri al, 
(59) 


where only the fermion part of the path integral has been 
written. This correlation function is non-zero provided 


|v1 — 24] ~ [aa = 24] S$ 1/m, (60) 


where m~! is the width of the ZM bound to the instanton 
at Z+. 


A nonzero value for such an anomalous (Gor’kov) 
Green’s function is suggestive of symmetry breaking. It 
is indeed a gauge-invariant object, since %ı and %2 cou- 
ple to the dynamical U (1) gauge field a, with opposite 
charge, but transforms nontrivially under the global U (1) 
symmetry associated with boson number conservation in 
the original microscopic model (see Sec. II). To investi- 
gate whether symmetry breaking indeed occurs, we add 
a weak symmetry-breaking source J and consider the 
fermion part of the path integral in Eq. (54), 


Zrla}, J] = J Dla e BE ea, (61) 


where the source J (x,y) is generically nonlocal with some 
spinor structure. This is not a valid term by itself, since 
it renders the Hamiltonian non-Hermitian (or the action 
non reflection positive). However, the Q= —1 sector will 
provide the needed conjugate term (Sec. V C). In nonlo- 
cal expressions like the source term in (61), Wilson lines 
should be inserted to maintain local gauge invariance. In 
accordance with our neglect of fermion-photon interac- 
tions at this stage, and because the final form of the ’t 
Hooft vertex will be a local interaction, we do not write 
down these Wilson lines explicitly. Working to linear or- 
der in J, the action can be simplified using the mode 
expansions in Eqs. (56)-(57) as 


Sela, J] = Sela] + f (ev) 
x >u (a — z})J(z, y)v;(y = z4 )niXj 
= Splat, J] + [tenezm 


£ f 
+) (uh Jv;)noxs + D5 (uf Joo), 
J i 


(62) 


where d3(x,y) = d?zx dy, and all the non-ZM contribu- 
tions have been subsumed into Splat, J]. Using the func- 
tional measure (58), 


Zp(af, J] = J dno dxo D' (7, n, X, xe Srei 
x fi- f a(x, v)(uf Fe0}nox0| 
x ji- 5 f a°(e.s)(u5.7e.)m 
x j- D I dæ, yu] Jv), (63) 
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Figure 1. Fermion pair annihilation due to the source J (x, y) 
and instanton at z+. 


the square brackets coalesce into the expression 


1- f v) 


/ / 
x (us Jvo)noxo— 5 (ud Jui)noxi— 5 (u Jvo)nxo] 


/ 
+> J EE DEE AI Joox: 
ij 
(64) 


All the terms except the second vanish because of the 
functional measure (58); dyo kills the third, dno the 
fourth, both kill the first, and dizo dXizo kills the last. 
Therefore, we obtain a nonvanishing result: 


Ze(at, J)= | edy ug (e-z) (e, y)volu-2s)K, (65) 


where K denotes the contribution of nonzero modes. The 
order of integration, dno dXo, has yielded a minus sign. 
The factor K is independent of J when working to first 
order in the weak source J. Recalling that the ZM spinors 
uo, vo have radial parts of the form exp(—m|a—z+|), sim- 
ilar to the free-fermion propagator, we note that Eq. (65) 
resembles the structure of the Feynman diagram in Fig. 1. 


As an ansatz for the result of integrating out instan- 
tons, we are thus motivated to consider the path integral, 


I[J] = J Dap T eaa 


x f čarda Algauw] FT (xr), (66) 


where Sp (written without a source argument) is the free 
Dirac action in the absence of instantons, A is some con- 
stant to be determined, and w1,2 are spinors (possibly 
spacetime dependent) to be determined. A in the in- 
sertion can pair up with a w in the source term to give 
the free propagator that we desire. A,w are malleable 
quantities that must be fixed to obtain the exact result 
in Eq. (65). The specific form of the insertion is moti- 
vated in hindsight by the calculation that follows. Taylor 


expanding the source exponential, we obtain: 
= f Dba, baler 


x |1— | d(x, yyy] (x) J re) 


f Ena) Atba(r2)wowt wy (21) 
= —Al(0 o) fa T, Y, £1, £2) 
x (YP @)I A (x, yb WPI eaog), 
(67) 


defining (O(Ya,Va)))= ri S PPa Vale" O(a, ba). 
The (only) effect of I[0] is to restrict Wick contractions 
to connected diagrams, and therefore 


x 


=A | èle, y 21,2) 
x [Gi (x — x1)wı]" J (x, y)[G2(y — z2)w2], (68) 


where 


G(x — y) = —(T4(z)¥(y)),» (69) 


is the free, Euclidean, Dirac propagator. 
Comparison of Eq. (68) with Eq. (65) tells us Z[J] = 
Zrla}, J] if we make the identifications 


A=K, 


[eacile — £1)w1 = —ug(x — z+), 
fë — z2)w2 = —vo(y — z+). (70) 


The minus signs on the right-hand side are conventional, 
and I[J] = Zr[a{,J] even without them. Clearly, the 
second and third equalities demand w; =w;(a;— 2+). A 
shift of integration variables x; > x,+7z reveals that 
these are Fredholm integral equations of the first kind, 
with solutions 


w= —Gy ‘uo = 


W2 = -G3 wo = 


(J + m)uo, 
(Ø — m)vo. (71) 


Substituting the results for A and w into Eq. (66) gives 
je Sef (ey) VT Ih fë ats 


J= [Dar 


x bo(x2) [Kw (22-24 )w] (x1 —z+)] PF (21). 
(72) 


The single instanton of the Q = 1 sector has been inte- 
grated out. This path integral still evaluates to zero if 
the source J is switched off, since the free Dirac action 
Sp conserves fermion number. So far, all we have done 
is rewrite zero in different garb. 
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C. Q=-!1 anti-instanton sector 


A similar analysis can be carried out for the Q = —1 
sector, assuming a single charge —1/2e instanton (anti- 
instanton) localized at z_ and described by a background 
gauge field a,,. The final result can actually be written 
down by the requirement of reflection positivity alone, 
but for the sake of completeness we explicitly rederive 
the result here. 

The anti-instanton sectors gift the adjoint Dirac oper- 
ators Dİ with ZMs, and thus contribute to anomalous 
Green’s functions of the form (paq). This means one 


has to add a source term of the form p2 JY] and consider 


a= f Dia, Vale 


where J = ozJto, with J being the source added in the 
discussion of the Q = 1 instanton sector. This is deter- 
mined by the requirement of reflection positivity of the 
action with the full source term (YI Jw2+y2JVT). 


The mode expansions are now 


—Sprlaz]—f a° (£y) Jp] , (73) 


wa) = D ule- 2) 


palz) = 5 vi(x — 2—)Xis 
p(x) = x)io 7 2 U;(z x)» 


p(x) = £)Xo + pe vi(z_—2)Xi, (7A) 


where u; = &] and ŭo is the ZM of DÝ calculated in 
Sec. IV C. Inserting this into Zp|a,, J] gives the analog 
of Eq. (65) for an anti-instanton as 


Zp(a,,, J] =-K | aPed*y io(2 2) Fle, yap (2-—v) 


=-K favea’y 0g (a—z_)J (x, y)tio(y—z_). 
(75) 


The Feynman-diagram interpretation of Fig. 1 holds 
again but for fermion pair creation due to an anti- 
instanton at spacetime coordinate z_. There is an extra 
minus sign here compared to Eq. (65), due to the or- 
der of the measure dno dxo. This is important to obtain 
a reflection positive action at the end. The contribution 
from nonzero modes, subsumed into K, is the same as the 
Q=+1 sector if the eigenfunctions of DİD and DD? are 
used to account for the non-ZM contributions in mode 
expansions of fermion fields. This is because the nonzero 
eigenmodes of both operators are paired with the same 
eigenvalues. In any case, the precise numerical factor is 
unimportant here. 


Similar to the analysis in the previous section, we con- 
sider with an insertion a path integral 


TJ] = f Datas a? (x,y)b2 IVT 
x f (es, 12)AUT (m)l), (76) 


which can simplified as before to 


TJ] = -A f èle, y 21,2) 
x [E]Ga(x2 — £)]J (æ, y) [G] (a1 — y)&J- (77) 


Demanding equality with Eq. (75) sets A= K and pro- 
vides equations to solve for £1 9: 


[®ecten — y) (xı — z_) = tio(y — z_), 
[ene (vq — x)€o(aq — 2) = ŭols — z—). (78) 


Again, these are Fredholm integral equations of the first 
kind, with solutions 


& = —(-9+ m)iio, 


The explicit forms of the ZMs, quoted here from Sec. IV, 
are 


2 v2m 1/2 


“~~ r E Y4/2,0,08: p) = vo, 
2 V2M mr 
E Wane p) = Uo, (80) 


where uo, vo are the respective ZMs of D1, D2. The equa- 
tions (79) for ¿&; thus become 


& = (J — m)vo = we, 
E2 = (g + m)uo = W1; (81) 


where w; were defined in Eq. (71). Therefore, the final 
result for the fermion path integral in the Q =-—1 sector 
is 


Zrla;, J] = f ate d° (wy) ba Jh] fE 
x yi (z1) [Kwz(z1—z-)w] (2—2_)] Y2 (x2). 
(82) 


D. Resummation and a local Lagrangian 


Inserting the results of Sec. V B and V C into the par- 
tition function (49) of the full theory, where only Q = +1 
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instantons are kept, we obtain: 


Z[J] = J P@artha)De e7 S Pa Lo- S a? (x,y) (pI JbztbeI$T) 


love) AK N N 
x 5 ( a TL | ta f zay 
N=0 k=1 
x [e Pr ENT (awaz zu) (Yzel) 


tetro lataa (ar)wo(a—24)u (uzt), 
(83) 





where 





Lo = Y (J — id + mJy + Y2(9 + id — m)y 4 a (dc), 
(84) 


is the Lagrangian of Eq. (46) but absent instanton ef- 
fects, with the Maxwell term dualized. We have also re- 
instated fermion-photon interactions to maintain explicit 
gauge invariance. The k-product in Eq. (83) just gives 
the Nt power of the insertion and, summing over N, an 
exponential is born. Exponentiating and then setting the 
source J=0 results in a nonlocal effective action 


Seg = fèrt — AK f azdad®y 
x [ersBrele) +02 (2)un(a—2)uF (y—2)ubal) 


+ etro l) tO) (x)wa(£z— z)wt (y—z)v] (y). 
(85) 


Can this action be approximated by a local one? Because 
the ZM wavefunctions decay exponentially in spacetime, 
it is reasonable to expect so (recall Fig. 1 and the discus- 
sion surrounding it). A change of integration variables, 
xz—>zxr+z and y—>y+z, allows the rewriting of one of the 
terms in the ’t Hooft vertex (i.e., the instanton-induced 
action) as: 





ASer = -AK / dP zd? rd?y e227 (2) +44) 
x YT (x + z)wo(x)w] (y)do(y + 2). (86) 


Since wı and w2 are proportional to the radial part 
(e-™"/r) of the ZM, the dominant contributions to the 
x and y integrals are from small neighborhoods of «=0 
and y=0. Taylor expanding the Fermi fields in powers 
of x and y to leading (zeroth) order gives 


ASert x% -AK f as e272 (z)+4(2)] 


<W ( f areun(a) f yuto) vale) 
(87) 


Using Eq. (71), and denoting vo(p) = 


f Erw) = f&a (A — m)vo(z) 


z = im (ip — m)vo(p) 


= on [ Paw volz 


> —Vims?? f E 0,0 


"aa ey (88) 
and similarly, 
J d3a wi (x) = I d’z( + m)uo(2) 

= Vam f dr re™™ "faf e T 


"o ie (89) 


we find that the quantity appearing in brackets between 
pI (z) and w2(z) in Eq. (87) is 


f a e~P* up (z), 








[ae w(x) I d’ywT (y) œx -4 (e, + icy), (90) 


where the proportionality constant is some number which 
shall be subsumed into K. This implies the effective ac- 
tion is specified by the local Lagrangian 


e204 8/2) NT (o, +ioy) 2 


2 462 
Ke-™ /4¢ 
E aed 


+ e@milo+9/2) 5 (a, +ioy) Pt |. (91) 


Using the transformation © in Appendix B, with addi- 
tionally O(o(a)) = o (0x) for a scalar field, one can check 
explicitly that the local ’t Hooft vertex thus derived is 
reflection positive, and thus corresponds to an interac- 
tion that preserves unitarity of the underlying real-time 
quantum field theory. 


VI. PARTONS AND SYMMETRY BREAKING 


The original parton gauge theory had 2N;y=4 fermion 
flavors, described by the Lagrangian (6). In the preced- 
ing Sec. V, to simplify the calculation of the ’t Hooft 
vertex, we retained only two fermion flavors while pre- 
serving the U(1) global and gauge symmetries. It can 
be seen from the calculations in that section that instan- 
tons in CQED3 with 2Ny flavors of fermions with the 
given mass and charge assignments will induce a ’t Hooft 
vertex with 2N,y fermion operators. For example, in the 
case of the original Lagrangian (6) with four fermion fla- 
vors {w14,W2+4}, one considers the four Euclidean Dirac 
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operators: 


Dig = V0, — iy" ay, +m, 
Dot = APD 








+ iy" ay, — m. (92) 


This simply doubles the number of fermion ZMs in each 
instanton charge sector. For instance, in the Q=1 sector, 
each of the four operators possesses a ZM, and the mode 
expansions of the fermion fields now become: 


Wit (x) = uols i z4)Not Sl > u(x = Z4)Mik; 


= vo(z — 24)xor + D> v(x — z+)Xiz. (93) 


i 














Wot (x) 


This results in a path integral measure: 





= dng dng dxg dxo 
x D' (f+, n+)D'(X+,x+), (94) 


where the primed measure includes contributions from 
nonzero modes in the expansion (93). Therefore, to ob- 
tain a nonzero path integral, a four-fermion insertion of 
the form of], y2+¥]_2_ is required. Repeating the cal- 
culation of the ’t Hooft vertex in Sec. V with a four- 
fermion source and insertion yields an instanton-induced 
term (the exponentiated insertion) in the Lagrangian of 
the form: 


eM ei yT bos pI yp- + H.c., (95) 


D(di4,dr14, Pot, Wot) 



































where y=o, +io,, and we have used w=w'o, to rewrite 
the second term as the Hermitian conjugate of the first. 
Because w+, and 7_ create excitations with lattice mo- 
menta near the Dirac points K} and K_, respectively, 
the presence of an equal number of y, and w_ fields in 
(95) guarantees the ’t Hooft vertex respects the micro- 
scopic translation symmetry (since K4} + K_=0 modulo 
a reciprocal lattice vector). 

In the absence of instantons, (noncompact) QED; 
has a global topological U(1)iop symmetry associated 
with the conservation of the topological current jp? = = 


7A —€uvava. In the dual formulation, this symmetry is a 
shift symmetry of the dual photon o —>o+a, manifest in 
the Lagrangian (84). The parton theory with noncom- 
pact gauge fluctuations thus has the global symmetry 
U(1) x U(1)top, where the first U(1) is the boson num- 
ber conservation symmetry under which w+ > e?y14 
and 2+ — Y2+ [recall the choice of global charge as- 
signments in Eq. (6)|. The ’t Hooft vertex (95) shows 
that instantons have the effect of explicitly breaking this 
U(1) x U(1)top symmetry to a diagonal U(1) subgroup 
under which 























eee 
T 


pi> e pit, Yor tre, (mod 1). 


(96) 

















The latter transformation makes clear the fact that ø is 
a compact scalar field of compactification radius 1. The 


diagonal U(1) symmetry (96) is to be understood as the 
correct incarnation of the unique microscopic U(1) boson 
number conservation symmetry in the low-energy par- 
ton theory with compact gauge fluctuations, i.e., where 
Polyakov instantons are accounted for. 

Although instantons have been explicitly taken into 
account in the derivation of the four-fermion ’t Hooft 
vertex (95), the resulting effective theory is still an inter- 
acting gauge theory, and its infrared fate not altogether 
obvious. A natural route to confinement—our primary 
focus—is the instanton proliferation scenario, whereby 
the coefficient of the *t Hooft vertex (95) is assumed to 
run to strong coupling under renormalization group flow. 
One then expects spontaneous breaking of the global 
U(1) symmetry (96), with o acquiring an expectation 
value [38, 43, 59]. The ø field itself is the Goldstone mode 
of the broken continuous symmetry, and the microscopic 
hard-core boson system becomes superfluid, as discussed 
in Ref. [44]. Additionally, (95) shows that a constant 
0 parameter can be given a natural interpretation as a 
global shift in the phase of the condensate. 

However, (7) £0 only implies that the U(1) symmetry 
is broken to a Zə subgroup under which i+ > —Y1+ 
(8=7), since the ’t Hooft vertex contains two p, fields. 
In terms of the original constituent bosons, this cor- 
responds to a boson pair condensate (b(a)b(a’)) 4 0 
without single-particle condensation, (b(a)) = 0, which 
preserves an Ising symmetry b(æ) + —b(a) (see, e.g., 
Refs. [60, 61]). In terms of the fermionic partons, the 
order parameter (¢)],yW24V]_72_) £0 is analogous to 
that for charge-4e superconductivity [62], but without 
concomitant Higgsing of the U(1) gauge symmetry since 
(95) is manifestly gauge invariant (recall that ~ and 72 
carry opposite gauge charge under the dynamical gauge 
field). 

The residual global Zə symmetry in such a paired su- 
perfluid can be further broken [60, 61], yielding a con- 
ventional superfluid phase with single-particle conden- 
sate (b(a)) #0. In the current context, this occurs if a 
gauge-invariant fermion bilinear condenses, (wWwW2) Æ 0. 
The various possible spinor/valley index structures of 
such a bilinear (suppressed here) allow in principle for 
both translationally invariant condensates and spatially 
modulated ones, i.e., supersolid phases. 








VII. CONCLUSION 


In summary, we have presented a nonperturbative 
study of monopole-instanton effects in a (2+1)D par- 
ton gauge theory featuring Dirac fermions coupled to a 
compact U(1) gauge field—CQED3. This parton gauge 
theory is meant to encapsulate the universal low-energy 
physics of hard-core lattice bosons in the vicinity of a 
multicritical point separating fractionalized phases, such 
as boson fractional quantum Hall states, and conven- 
tional ones. While the compactness of the gauge field be- 
comes irrelevant in fractionalized phases, which support 
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deconfined excitations, we focused on developing an ex- 
plicit understanding of the instanton dynamics that leads 
to confinement in conventional phases. As our first main 
result, we showed that in contrast to CQED3 with mass- 
less fermions—an effective gauge theory describing the 
Dirac spin liquid—CQED3 with massive fermions sup- 
ports Euclidean fermion zero modes exponentially local- 
ized on instantons. The localization length of the zero 
mode “wavefunction” is found to be inversely propor- 
tional to the fermion mass, which in hindsight elucidates 
the absence, first observed by Marston, of normalizable 
zero modes in massless CQED3. While we did not prove a 
rigorous index theorem guaranteeing the topological sta- 
bility of such Euclidean zero modes, they were found to 
be in one-to-one correspondence with Hamiltonian quasi- 
zero modes occurring in the context of monopole operator 
dressing in conformal field theories associated with spin 
ordering transitions out of the Dirac spin liquid. In such 
theories, a nonzero fermion mass arises when the theory 
is canonically quantized on the sphere, and the result- 
ing Hamiltonian quasi-zero modes can be understood as 
“massive deformations” of true zero modes protected by 
the Atiyah-Singer index theorem. 


As our second main result, we combined semiclassi- 
cal methods with our zero mode solutions to show by 
explicit derivation that instantons mediate an effective 
four-fermion interaction in the gauge theory, known as 
the ’t Hooft vertex. This effective interaction explicitly 
breaks a spurious U(1) x U(1)top symmetry of the clas- 
sical parton Lagrangian to a diagonal U(1) subgroup, 
corresponding to the physical boson number conserva- 
tion symmetry of the microscopic model. Under the 
further assumption of confinement via instanton prolif- 
eration, we found that the ’t Hooft vertex could natu- 
rally lead to two distinct superfluid phases: an ordinary 
single-particle condensate, but also a boson pair conden- 
sate without single-particle condensation, in which the 
global U(1) symmetry is only broken to Zə. 


Looking ahead, our approach based on semiclassi- 
cal instanton techniques could be used to complement 
the Hamiltonian monopole-operator dressing approach 
to confinement transitions out of the Dirac spin liq- 
uid [10, 11]. Song et al. rely solely on microscopic sym- 
metries and write down deformations of the conformal 
QED3 Lagrangian consisting of (dressed) monopole op- 
erator/fermion composites allowed by those symmetries. 
Alternatively, ’t Hooft vertices containing similar physics 
could be explicitly derived as follows. In the two-step 
route to confinement advocated by Song et al. and men- 
tioned earlier, a fermion mass bilinear acquires an ex- 
pectation value before instanton proliferation proceeds. 
Applying the semiclassical methods employed here af- 
ter the first step, Euclidean zero modes for the resulting 
massive fermions could be searched for and used to de- 
rive a ’t Hooft vertex that would encapsulate the range 
of symmetry-breaking phases made possible by instan- 
ton proliferation. Finally, the proof of an index theo- 
rem for massive Dirac fermions in 3D Abelian instanton 


backgrounds would be a desirable extension of the results 
presented here. 
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Appendix A: Monopole miscellanea 
1. Monopole harmonics 


This appendix collates some well-known results on the 
theory of Dirac monopoles and is mostly self-contained. 
We first consider a spinless charge e in the field of a static 
point monopole at the origin, 

B= ah (A1) 
described by a Wu-Yang vector potential A [see Eq. (27)]. 
Classically, the spherical symmetry of the problem sug- 
gests conservation of angular momentum. The natural 
guess l=rx(p—eA) = rxmv, by minimal coupling, does 
not work because 


dl 
in at x mv), 
=r X mi, 
=r xe(v x B), 
eg 


=rx 73 (UXT); (A2) 
=l xr. (A3) 
mr 


This is generically non-zero, suggesting lis not conserved. 
Using a formula for the vector triple product, Eq. (A2) 
can be rewritten as [63] 


dl eg ,. 
Gots wat xr), 
eg : 
73 [re —(r- t)r] f 
ee | ap 10r 
~ i j 2 dt |. 
1 1 
= eg zr — ltr , 
d r 
= g he7) Aa] 
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This implies a conserved angular momentum 


L =r x (p-—eA)-— qf, q = eg € Z/2. (A5) 
One can explicitly prove (post-quantization) that 
[L;, L;] = l€ijkLk- (A6) 


Since [r?, L] = 0, these two operators can be simulta- 
neously diagonalized and L can be studied for fixed r. 
Also, since [L,, L?] = 0 and [L;, Lj] = icijkLk, we have 
the familiar 


LY, pil 
L2Yq,1,u( 


+9) = L(L + 1)¥q,r,u (0, p), 


0 
0, p) = MY, L,m (0, p). (A7) 
The sections Y} ,m are called monopole harmonics, and 
their exact form is gauge dependent, which means north- 
ern and southern versions differ by a gauge transforma- 
tion in the Wu-Yang formulation. Only L, M are quan- 
tum numbers, while q is a parameter that determines one 
complete set of harmonics. Just based on the su(2) alge- 
bra, allowed values of L must be a subset of {0,1/2,1,...}, 


while Me{—L,—-L+1,..., L}. However, 


L? = (l— qt)? 
=P+q-o(l-t+2-0 
=P +g. (A8) 





For fixed 2q EZ, this gives a bound on the eigenvalues 


L(L+1) 2 &@. (A9) 


The solution of the inequality above is L > |q|. To prove 
this, we may take q > 0 without loss of generality as 
the inequality is independent of sgn q. Factorization and 
substitution of q=n/2, where nEN, gives 


(.- (z- EET) si i 


Both brackets must be of the same sign. Positivity of L 
implies 


s —1+v1 +n? 


L 
2 


L>n/2=q satisfies this inequality, since for n>0, 


—l1+vy1 +n? Z —-l+vVl+n?+2n n 
2 5 2 EP 





To see that this is the smallest satisfactory half-integral 
L, note that the next smallest value of L does not satisfy: 


n=1 _ —1 + vn? z -1l+vVJ1-4n? 
Q 2 = 2 f 





We thus have the result that L> |q|. 


Written in spherical coordinates, Eq. (A5) reads 





LN/S = —gt+6[i csc 00,.—q(cot OF csc 0)]— iða. (A11) 





This implies 





L, = -i0, F 4q, (A12) 
for the 2 component, which has northern and south- 
ern eigenfunctions of the form exp[i(M+q)y]. The re- 
quirement of a single-valued wavefunction then mandates 
(Mq) €Z, which is satisfied if M is (half-)integral when- 
ever q is (half-)integral. Together with L > |q|, this de- 
termines the allowed values of (L, M) as 








Le {lq la| +1,--}, Me{-L,-L+1,..., L}. (A13) 


For completeness, we provide a general formula for the 
monopole harmonics Y} zm in terms of the Wigner D- 
matrix. An elegant derivation can be found in Ref. [64]. 
In the northern hemisphere, 


2L+1 * 
Yq,t,.M(9n,%) = 4/ -i [Dir (20, -~)], (A14) 


where Oyn € [0,7). The southern versions (valid on the 
south pole) are obtained by a gauge transformation, 


¥,2,m (05, p) =e ?Y L, mM (Ons): (A15) 


The Wigner D-matrix is defined in terms of Euler angles 
(a, 8, y) as 


Diy mla, B, y) = (jm | e7 eie jm), 


ae T (A16) 


Using the formula above, the first two q=1/2 harmonics 
are given by 





1 o 

Yapan 9) = -o 
1 Ø 

(On, p) = -= cos 5 


V 2T 2” 


in the north. Their southern versions are given by a 
gauge transformation exp(—iy). 


Yı (A17) 


1i _1 
233233 9 


For g=—1/2, the first two northern harmonics are 


0 
(0N, p) = F= COS 5; 
22> T 2 


N= 

N= 

n= 
N 





(A18) 


; 6 
e`’? sin -, 
2 


with their southern versions now obtained by a gauge 
transformation exp(iy). 
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2. Monopole spinor harmonics 


We now consider a spin-1/2 particle of charge e in a 
monopole background. The total angular momentum for 
a spin-1/2 is 


1 
J=L+50. (A19) 


The allowed eigenvalues of J? and Jz, respectively de- 
noted j(j+1) and m,, follow from rules for addition of 
angular momenta: 


; 1 1 
je {L~1/2,b41/2}= fld-3 la+ 3} 
m; E€ {-j,-j + 1,...,j} (A20) 


The same rules also provide the (angular) eigensections 
Dms called monopole spinor harmonics, 


J+miYa jg a 
$ 


; 1 
Ve E | VA 
GJ,™5 V2) Vi- mY, j-2 mth 





yl? (9 p) = bf VIN MG Yai m 
Bpmg M242 | Sot MAG i ma J 
(A21) 


where Yj, are the monopole harmonics defined in 
Sec. A 2, and their coefficients are Clebsch-Gordan. For a 
given q, these spinor harmonics are a complete, orthonor- 
mal set of 2-spinor eigensections of L?, o?, J?, Jz. 

For use in the main text, we also record here the action 
of o-fon vë T which can be explicitly evaluated using 
Eq. (A21). Alternatively, since ø - ê commutes with J, 
the most it can do is mix the L= j+1/2 states. A general 
formula is 


(o s f) yit! = ayit? + by? 


qj, Mj E7 q,j,Mj qj,mj' 








(A22) 








Substituting in Eq. (A21) and using the known forms of 
the monopole harmonics provides linear equations for the 
coefficients, which turn out to be [65] 


2q 
2j +1 
(2j + 1)? — 4q? 
2j +1 l 





(A23) 


Appendix B: Self-adjoint operators 


This appendix elaborates on some technical aspects of 
the path integrals studied in Sec. V B-V C, particularly 
the difficulties involved in suitably defining the functional 
measure and connections with index theorems. We shall 
first proceed along a standard route used in the physics 


literature to define path integral measures. The path 
integral of interest is 
Zelof]= f Dyes to, (BI) 


where ay describes a single charge 1/2e instanton located 
at z}. Although an even number of fermion flavors are 
required for this theory to make physical sense, a sin- 
gle flavor is sufficient to highlight some of the general 
mathematical difficulties that arise in this problem. The 


massive) Euclidean Dirac operator D is defined as 
P 
D = "ô, — iy"ay +m, 


sied a, . (7 r 





3 
) z tm. (B2) 
7 
The (naively taken) adjoint of the Dirac operator is 
Di = =q" ð, + iqa +m 


Saoi a, . (7 r 





3) to . J +m. 
(B3) 

The second lines of Eqs. (B2)-(B3) are the results of 

Sec. IV B-IV C. These operators are not Hermitian, but 


one can consider the Hermitian combinations DİD and 
DD! with eigenvalue equations 


D'Du,(a — z4) = aju; (£ — 24 
DD ai, (x — 2.) = G,ti,;(a — z4), 


Ww 








(B4) 


for an instanton located at z+. 

We have used the term “Hermitian operator” to refer 
to what is called a “symmetric operator” in the mathe- 
matical literature on unbounded operators [66, 67]. Act- 
ing on a Hilbert space, a densely defined Hermitian 
(or symmetric) operator A satisfies (u, Av) = (Au, v), 
for any u,v € Dom(A). To be self-adjoint, that is for 
A = At to hold as an operator equation, one also re- 
quires Dom(A) = Dom(A‘), which does not follow from 
the Hermiticity condition for unbounded operators (such 
as the Dirac operator under consideration), and usually 
Dom(A') C Dom(A). Only self-adjoint operators have 
the desirable properties of possessing a complete set of 
eigenfunctions and real eigenvalues. However, this fact 
is typically ignored, and one proceeds to use the eigen- 
functions of DİD and DD’ as a basis to facilitate a mode 
expansion of the Fermi fields in the path integral, assum- 
ing such operators are indeed self-adjoint. This is usually 
harmless, but not so in current circumstances as we shall 
momentarily show. In any case, a loose mathematical 
justification of this ignorance can be made by assuming 
that the Hermitian operators above possess a unique, or 
a family, of self-adjoint extensions. An arbitrary choice 
in this family will have the required property of possess- 
ing a complete basis of eigenfunctions to facilitate mode 
expansions. 

However, a self-adjoint extension of a Hermitian opera- 
tor A involves imposing boundary conditions on its eigen- 
functions which effectively shrink or enlarge Dom(A) and 
Dom(A‘) until they are equal. The eigenfunctions of the 
self-adjoint extension will form a complete basis only for 
the final domain Dom’(A). Assuming DİD or DD? have 
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been made self-adjoint, mode expansions of Fermi fields 
in terms of the eigenfunctions of these operators effec- 
tively assume that the space of fields being integrated 
over in the path integral is the same as Dom(D'D) or 
Dom(DD"). This typically does not warrant close analy- 
sis since one hopes (usually correctly) that all important 
physical effects are accounted for in this procedure. In 
the present case, as we show below, the ZM solution lies 
outside these domains and is thus missed if the eigen- 
functions of DİD or DD" are used for a mode expansion 
or in the definition of the functional measure. 

We start with the following paradox. The ZMs of D 
and Di were calculated in Sec. IVB and IVC respec- 
tively. For g=—1/2e, the operator Dt has a normalizable 
ZM 


v2m 1/2 


vo (r, 0, p) = OY 72,0,0(8: p), (B5) 


which implies 


(Wo D'Y) = 0. (B6) 
However, using Eq. (B2), one finds explicitly that Dio = 
2mwo , which seems to imply 


ee ~ 2 
(Diy 5) = mlir || #0. (B7) 
The implied consequence is that (by ; Dipo ) F# 
(Dug Wo ). Before a resolution of this is pointed out, 
we note that Dig =2myp implies then that DÌ Dýg =0 
so that dim ker DİD — dim ker DD‘ = 0. It is the latter 
that is typically calculated as index(D) in typical proofs 
of index theorems used in physics. This highlights the 
difficulty in producing an index theorem for the current 
scenario, and also in defining the functional measure of 
the path integral using eigenfunctions of D'D or DD’. 

The resolution of the paradox lies in a careful exami- 
nation of the domains of the operators Dİ and D, which 
turn out to be subspaces of square integrable functions. 
In the subspace of spinors w with fixed angular part 
Yit y), and for eg =1/2 (i.e., in the Q=1 instanton 
sector discussed in Sec. V B), the action becomes: 


S= fèr wip, +m), (B8) 


where one defines a “radial momentum operator” [68], 


pe =-5@-V40-8) =-i(0, +2). (B9) 


The adjoint Dirac operator Dt in Eq. (B3) has been 
naively derived from the form of D by essentially assum- 
ing pr is Hermitian on Dom(D). However, this Hermitic- 
ity condition is violated on the ZM, 


1 
polr) x =e ™, 
r 


(B10) 


of the operator (ip, +m) appearing in the action (B8), 
for (pro, Yo) Æ (Wo, Pro). This is the reason for the 
paradoxical equations (B6)-(B7). However, to restrict 
the path integral over w to a function space on which 
Pr is Hermitian is to exclude ZMs and their associated 
physics. To determine the space of fields (Y and w) that 
one should integrate over, we use the necessary condition 
that the Minkowski action must be real-valued. 

The reality of the Minkowski action in a unitary quan- 
tum field theory translates to Osterwalder-Schrader or 
reflection positivity of the corresponding Euclidean ac- 
tion, which is invariance O(S) = S under a form of com- 
plex conjugation followed by Euclidean time reversal [69]. 
This transformation acts on fermions as an involution 
of the Grassmann algebra [70], which for our particular 
choice of Dirac matrices can be chosen as 


O(da(x)) = of" Up (6z), (B11) 
O(ba(x)) = of? bp (0x), (B12) 
where 0 flips the sign of the time (z) coordi- 


nate. Additionally, © complex conjugates c-numbers 
and reverses the order of Grassmann variables, e.g., 
O(waveypy) = O(,)O(vg)O(Wa). Gauge fields trans- 
form as O(ao(#)) = —ao (0x) and O(a;(x)) = a;(Ax) [71]. 
One can show that under O, Eq. (B8) transforms as: 


Tr=CoO 


(8) =5+ / do (=r By) "=, (B13) 


where the boundary term follows from an integration 
by parts. However, for reflection positivity to hold, the 
boundary term is required to vanish. 

The upper limit of the boundary term vanishes if we 
require all fields to be square integrable. Say y~r~°x as 
r—0O, where x is a Grassmann number without r depen- 
dence. Then square integrability requires lim, r*?° to 
exist, which means 6 <3/2. Since square integrable func- 
tions are therefore at most as singular as r—%/?+*, where 
€>0, the lower limit of the boundary term is at most as 
singular as 


r2 


lim -aeara (B14) 


where ô>0. The existence of this limit requires 6+¢>1. 
The choice e = ô > 1/2 restricts both field integrations 
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